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Report  Title 

Direct  Numerical  Simulations  of  Very  Stable  Atmospheric  Boundary  Layers 

ABSTRACT 

The  objective  of  the  research  is  to  address  a  number  of  outstanding  issues  related  to  very  stable  atmospheric  boundary  layers,  including 
criteria  for  the  collapse  of  the  turbulence,  the  intennittent  behavior  and  the  strong  anisotropy  of  turbulence,  horizontal  dispersion  and 
meandering  of  plumes,  and  the  local  scaling  of  the  velocities  in  terms  of  predictable  or  measurable  parameters. 

The  approach  to  the  research  is  direct  numerical  simulation  of  stably  stratified  Ekman  boundary  layers.  Very  high  resolution, 
three-dimensional,  time-dependent  simulations  are  carried  out  on  computational  grids  of  the  order  of  1,000  grid  points  in  each  spatial 
direction, using  highly  accurate  spectral  and  compact  numerical  methods. 

The  results  of  these  simulations  permit  a  complete  statistical  description  of  velocities, Reynolds  stresses,  and  other  terms  in  the  energy 
budgets.  The  effect  of  stratification  on  the  length-scales  of  the  flow  is  investigated  using  horizontal  spectra  of  various  important  dynamical 
quantities.  The  detailed  infonnation  from  the  DNS  allows  the  use  of  conditional  sampling  to  compute  the  statistical  properties  of  turbulent 
patches,  as  well  as  to  identify  their 

sources  (internal  waves,  low-level  jets,  etc.).  Finally,  the  simulation  results  will  also  be  useful  as  a  guide  to  field  measurements,  and  to 
perform  a  priori  tests  of  turbulence  models. 
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Scientific  Progress 


Our  results  use  the  technique  of  direct  numerical  simulation  to  study  stable  atmospheric  boundary  layers.  In  particular  we  have 
examed  when  and  how  the  turbulence  in  the  layer  essentially  shuts  down  due  to  night-time  cooling  at  the  ground,  as  observed 
in  the  CASES99  field  experiment.  Of  course  turbulence  in  the  atmospheric  boundary  layer  has  many  causes,  including  the 
aforementioned  ground-generated  turbulence,  as  well  as  turbulence  due  to  breaking  internal  waves,  turbulence  resulting  from 
the  shear  in  low  level  jets,  turbulence  due  to  topographical  features,  and  other  mechanisms. 

Direct  numerical  simulation  is  a  computational  technique  that  has  been  used  increasingly  in  the  study  of  turbulence  flows,  and 
now  often  in  the  study  of  turbulent,  stably-stratified  flows.  In  direct  numerical  simulations,  all  of  the  relevant  length  and  time 
scales  of  motion  are  resolved.  Because  of  the  limited  resolution  of  computer  simulations,  however,  in  this  approach  only  flows 
with  Reynolds  numbers  well  below  those  in  the  atmospheric  boundary  layer  can  be  treated.  Therefore,  only  problems  and 
issues  that  are  not  strongly  Reynolds-number  dependent  can  be  addressed  by  this  approach.  Since  the  'local'  Reynolds 
number  is  relatively  small  in  the  lower  part  of  the  surface  sublayer  under  strong  stratification,  the  collapse  of  ground-generated 
turbulence  is  exceptionally  well  suited  for  analysis  with  this  technique. 

The  results  of  our  simulations  using  smooth,  flat  ground  suggest  that  ground-generated  turbulence  shuts  down  when  the 
Monin-Oboukov  scale  reduces  to  the  order  of  the  buffer  region  in  the  boundary  layer,  which  region  is  located  very  close  to  the 
ground.  Results  from  non-stratified  boundary  layer  studies  have  indicated  that,  over  flat  surfaces,  the  buffer  region  is  the  home 
of  a  nonlinear  cycle  that  sustains  the  turbulence.  Our  results  indicate  that,  once  the  stratification  is  able  to  sufficiently  disturb 
that  cycle,  ground  generation  of  turbulence  is  shut  down.  When  these  arguments  are  extended  to  more  general  rough 
surfaces,  the  suggestion  is  that  the  turbulence  generation  over  rough  surfaces  shuts  down  when  the  ratio  of  the 
Monin-Oboukov  length  scale  and  the  roughness  sublayer  height  become  of  order  one.  This  suggestion  is  consistent  with  the 
results  from  CASES99,  and  also  suggests  criteria  to  be  used  in  RANS  and  large-eddy  simulations  of  atmospheric  boundary 
layers.  A  paper  based  upon  our  results  has  been  published  in  Boundary  Layer  Meteorology.  Also,  direct  numerical  simulation 
results  were  supplied  to  Dr.  Keith  Wilson,  ERDC-CRREL  for  use  in  testing  his  modeling. 
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